Draft version December 11, 2012 

Preprint typeset using I4TgX style emulateapj v. 8/13/10 



BROKEN AND UNBROKEN: THE MILKY WAY AND M31 STELLAR HALOES 
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ABSTRACT 

We use the iBullock & JohnstonI (I2005h suite of simulations to study the density profiles of L*-type galaxy 
stellar haloes. Observations of the Milky Way and M31 stellar haloes show contrasting results: the Milky Way 
has a 'broken' profile, where the density falls off more rapidly beyond ^ 25 kpc, while M3 1 has a smooth 
profile out to 100 kpc with no obvious break. Simulated stellar haloes, built solely by the accretion of dwarf 
galaxies, also exhibit this behavior: some haloes have breaks, while others don't. The presence or absence of 
a break in the stellar halo profile can be related to the accretion history of the galaxy. We find that a break 
radius is strongly related to the build up of stars at apocentres. We relate these findings to observations, and 
find that the 'break' in the Milky Way density profile is likely associated with a relatively early (~ 6 — 9 Gyr 
ago) and massive accretion event. In contrast, the absence of a break in the M3 1 stellar halo profile suggests 
that its accreted satellites have a wide range of apocentres. Hence, it is likely that M3 1 has had a much more 
prolonged accretion history than the Milky Way. 



1. INTRODUCTION 

A diffuse envelope of stars surrounds the Milky Way 
Galaxy. This halo of stars only contributes a meager few 
percent to the total light, but comprises the oldest, and most 
metal-poor stars in the Galaxy. The dynamical timescales in 
the radial range of the halo stars 10 — 100 kpc) are long, 
so these stars can preserve their initial conditions. It is widely 
recognized that by studying the phase-space and chemical 
properties of the stellar halo we have the opportunity to un- 
ravel the accretion history of our Galaxy. 

The density profile of the Milky Way stellar halo has been 
studied extensively over the past decade. Early work extend- 
ing out to 25 — 30 kpc found that the halo follows an oblate, 
single power-law distribution with minor- to-major axis ra- 
tio q ~ 0.5 — . 8, and power-law index a ^ ^ ~ 4 (e.g . 
lYanny et al.l2000l:lNewberg & Yannyl2006tlJuric et al.l2008li . 
More recent work probing further out in the stellar halo has 
found that the stellar density falls-off more rapidly beyond 
~ 25 — 30 kpc with a pow er-law index of ann t ^4 — 5 be- 
yond the break radius (e. g. IWatkins et al.ll2009l iDeason et al.l 
1201 Ibt ISesaret al.ll20IIh . One could argue that the density 
dist ribution can just as eas ily be described by an Einasto pro- 
file (lEinasto & HaudI 19891) . which allows for a steeper fall-off 
at larger radii, without resorting to a break. However, whether 
a broken power-law or Einasto profile are preferred, in both 
cases there exists a characteristic scale (break radius or effec- 
tive radius) which deserves a physical justification. Is a break 
radius a ubiquitous feature of stellar haloes? Unfortunately, 
the low surface brightness of stellar haloes inhibits the study 
of individual galaxies beyond the local group. However, our 
nearest neighbor, M31, provides a useful comparison. 

Early studies called into question whether a stellar halo 
even exists in M31: out to _R ~ 30 kpc the stellar density 
profile seems to be a continuation of the M31 bulge (e.g. 
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iPritchet & van den Berghl 119941: iDurrell et al.l 12004 . How- 
ever, more recent work extending to larger projected dis- 
tances (out to R ^ 100 — 200 kpc) does indeed find evidence 
for a metal-poor halo compon ent with a power-law index of 
a ^ 3.3 (e.g. iGuha thakurta et al.l I2005|: [ Irwin et all 



llbata et al.l 12007 



^urteau et all 1201 It [Gilbert et all 
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However, there is no evidence for a break in the M3 1 stellar 
halo profile: the density profile follows a continuous single 
power-law all the way from i? 30 kpc to at least i? ^ 90 
kpc. 

Of course, any measure of the smooth, underlying stellar 
halo density profile (if indeed, it does exist), is hampered by 
the presence of substructure. A wealth of subst ructure has 
now b een discovered in bo th the Milky Way (e.g. llbata et aTl 
1995; Newberg et al. 2002; Be lokurov et al. 2006; Jur ic etalj 
2008; Belokurov et al. 2007) and M3I ( e glbata et al. ll2001l: 
Ferguson et al. 2002; Gilbert et al.ll2007l: iMcC onnachi e et alj 
2009i) . These discoveries have strongly affirmed the theoret- 



ical predictions from numerical simulations that the majority 
of the stellar halo is built up from the accretion products of 
satellite galaxies. 

Numerical simulations of stellar hal oes have developed 
substa ntially over the past few years. IBullock & JohnstonI 
(120051) . hereafter BJ05, presented a suite of 1 1 stellar haloes 
built enti rely from the disrupted of accreted satellites. More 
recently, ICooper et al.l (120 1 Ol) used the Aquarius simulations 
and a dark matter particle 'tagging' metho d to produ c e stel- 
lar haloes in a fully cosmological context. IB ell et all (120081) 
studied the 'lumpiness' of main sequence turn-off (MSTO) 
stars in the Milky Way stell ar halo, and found consistency 
with the IBullock & JohnstonI (12005) simulations. This led the 
authors to suggest that the Milky Way stellar halo is consistent 
with b eing built up entirely by ac cretion. However, Xue et al] 
( I2OI II) and lDeason et al.l ( 1201 Ibl) find a somewhat 'smoother' 
stellar halo when traced by blue horizontal branch (BHB) 
stars. These observations perhaps suggest accretion may not 
be the only formation mechanism of the stellar halo. Recent 
hydrodynamic cosmological simulations postulate that some 
fraction of t he inner stellar halo is made up of s tars formed 
in situ (e.g. IZolotov et al. 20091: iFont et akfeoill) . However, 
stars can diffuse relatively quickly in configuration space (as 
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opposed to velocity space), so the spatial structure of the stel- 
lar halo can be smooth, even if it is built up from merging and 
accretion. Thus, a relatively smooth stellar halo could also 
signify an early accretion history. 

In this paper we study the density profiles of stellar haloes 
drawn from the BJ05 suite of simulations. We investigate the 
origin of 'broken' density profiles, and in particular, address 
why some haloes have an obvious break (e.g. the Milky Way) 
and why some do not (e.g. M31). We know that accretion is 
at least an important (if not the sole) contributor to the stellar 
halo. By comparison with the BJ05 simulations we hope to 
link the observed stellar halo profiles of local galaxies to their 
accretion histories. 

The paper is arranged as follows. In Section 2 we briefly 
describe the BJ05 simulations. In Section 3 we study the stel- 
lar halo density profiles of the 1 1 haloes drawn from the BJ05 
simulations. In Section 4 we investigate the origin of broken 
halo profiles and we discuss the implications for the Milky 
Way and M31 stellar haloes in Section 5. Finally, we summa- 
rize our main findings in Section 6. 

2. BULLOCK AND JOHNSTON SIMULATIONS 

The Bullock & Johnston simulations (see BJ05 for full de- 
tails), are a suite of 11 high resolution stellar haloes built 
up from the accretion of dwarf galaxies onto a Milky Way- 
like potential. The parent galaxy is represented by a time- 
dependent, analytical potential consisting of bulge, disk and 
halo components. The accretion history of each halo is ran- 
domly drawn within the context of a ACDM Universe; each 
accretion event was assigned a binding energy and orbital 
eccentricity drawn from the orbital distributions of satellites 
observed in cosmological simulations of structure formation. 
High resolution N-body simulations were then run to track the 
evolution of dark matter particles in each accretion event. At 
the present day (z = 0), individual haloes are built up from 
the disruption of > 100 dwarf galaxies. 

The baryonic component of each dwarf galaxy is followed 
by assuming that the cold gas inflow tracks the dark matter 
accretion rateQ. The gas mass is then used to determine the 
instantaneous star formation rate and to track the buildup of 
stars within each halo. The stellar matter in each dwarf galaxy 
is associated with the most tightly bound material in the halo. 
Variable mass-to-light {M /L) ratios for each particle are as- 
signed in order to produce realistic dwarf galaxy stellar pro- 
files; the M/L is chosen so that the luminous matter in the 
infalling satellites initially follows a King model embedded 
within an Navarro-Frenk- White (NEW) dark matter potential. 
The King profiles were chosen to specifically reproduce struc- 
ture (e.g. sigma vs. core relations) observed for dwarf satel- 
lites today. The time dependent chemical con tent of the dwarf 
galax ies is also included in the models (see [Robertson et al.l 
12:0051 andlfont et al. 2006 for more details). The stellar pop- 
ulation models include enrichment from Type la and Type II 
supernova, as well as stellar winds. A physically motivated 
supernova feedback prescription is used to reproduce the lo- 
cal dwarf galaxy stellar mass-metallicity relation. The rate of 
chemical enrichment is calculated analytically, and the abun- 
dances of H, He, Fe, and the a-elements O and Mg are tracked 
in the simulations. 

BJ05 verified that their models are able to reproduce a 
number of observational constraints, such as the number and 

^ This assumption holds prior to infall into the parent halo; the star forma- 
tion is truncated soon after each satellite is accreted. 
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Table 1 

The global properties of the 1 1 BJ05 stellar haloes. We give the 
halo ID from galaxia (see http://galaxia.sourceforge.net/) and BJ05 
in parenthesis, the best-fit single power-law slope, the best-fit 
broken power-law slopes and break radii and the KS test (SPL 
versus BPL) probability. 



distribution of structural parameters of the Milky Way's satel- 
lite population. Furthermore, the resulting stellar haloes have 
a similar luminosity (^ lO^L©) and density profile to the 
Milky W ay {a ^ 3.5), an d contain a similar amount of sub- 
structure jBell et al.ll2008l) . The scatter between the 1 1 simu- 
lated stellar haloes is due to their different accretion histories; 
we will explore the effect of accretion history on the stellar 
halo density profiles in the following sections. 

3. STELLAR HALO DENSITY PROFILE 

We investigate the density profile of halo stars belonging 
to the eleven parent haloes discussed in BJ05. Table 1 in 
BJ05 shows that the majority (80-90 %) of the stellar halo 
is built up from the 15 most massive accreted satellites. For 
this reason, we only consider the contribution to the stellar 
halo from the 15 largest satellites. We fit single power-law 
(SPL) and broken power-law (BPL) density profiles to stel- 
lar haloes between 10-100 kpc. To fit the density profiles, we 
com pute the density in radial bins and use the MPFIT IDL rou- 
tine ( lMarkwardtll2009l) to find the minimum Chi-square SPL 
or BPL density profile. The best-fit parameters for the density 
profiles are given in Table [T] In addition, we also fit profiles 
for stars belonging to individual accretion events. Provided 
the accreti on event is not too rec ent (less than 4.5 Gyr ago, cf. 
Figure 3 in lJohnston et alj|2008h . a BPL is a good description 
of the density profile (see e.g. Fig. [3]). 

The density profile for each stellar halo is shown by the 
thick black lines in Fig. [T] The dashed blue line shows the 
best-fit BPL profile for the stellar halo, and the break radius 
for this fit is shown by the vertical dotted line. The thinner 
lines show the profiles for the 5 most massive accreted satel- 
lites. The fraction (by mass) contributed to the stellar halo 
(between 10-100 kpc) decreases from thicker black lines to 
thinner red lines. 

Fig. [T]shows that a broken power-law profile is generally a 
good description of the stellar halo density profiles, where the 
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Figure 1. Radial density profiles of the 1 1 BJ05 stellar haloes. The thick black line shows the overall profile for the 15 most massive satellites 
that contribute to the stellar halo within 100 kpc. The thinner lines show the profiles for the 5 most massive accreted subhaloes. The fraction 
(by mass) contributed to the stellar halo decreases from thicker black to thinner red lines. The blue dashed line shows the best-fit BPL profile 
for the stellar halo, and the vertical dotted line indicates the corresponding break-radius. 



stellar density falls off more rapidly beyond the break radius. 
S uch 'broken' s t ellar h aloes are also seen in the stellar haloes 
of ICooper et al.l ( 120101) (see also BJ05). However, it is worth 
noting that some stellar haloes can just as easily be described 
by a single power-law. 

4. ORIGIN OF BROKEN PROFILES 

The stellar halo in the BJ05 models is a superposition of 
stars stripped from several satellite galaxies. Thus, to under- 
stand the origin of the global stellar halo profile, we must first 
investigate the profiles of stars belonging to individual accre- 
tion events. The density profile of stars stripped from an in- 



dividual satellite galaxy is well-described by a BPL (see e.g. 
Fig. |3]l. What causes this BPL profile? 

For each star particle, we find the apocentre and pericen- 
tre of its orbit. Assuming a spherically symmetric, stationary 
potential, we can estimate t he apocentre and pericentr e from 
the following equation (see iBinney & TremainSi 19871 Chap- 
ter 3): 

. 2[<l>il/u)-E] _^ 



= 0. 



Here, u ~ 1/r and the roots of this equation give the apocen- 
tre and pericentre. This equation can be solved for each star 
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particle given the potential and particle properties defined at 
redshift z = 0. This is a good approximation to the orbital 
properties of the stars at the time of stripping. However, the 
orbital properties may have evolved since the time of accre- 
tion onto the parent halo. In Appendix|A] we verify some of 
our deductions made from the halo star properties at redshift 
z = by tracing back the orbital histories of the accreted 
satellites. 

From the inferred orbital properties of the stars, we can find 
the average apocentre and pericentre of the stars belonging to 
an individual satellite accretion event. As the star particles 
can have different masses (see BJ05 for details), we compute 
the mass weighted average. In the top panel of Fig. [21 we 
show the average apocentre of star particles belonging to one 
satellite against the best-fit break radius for the density profile 
of these stars. The satellites from all eleven parent haloes 
are shown in this plot. We also show the average properties 
for the global stellar haloes with the star symbols. There is 
a clear correlation between apocentre and break radius; here, 
the dotted line indicates a one-to-one relation. 

The colors in the top panel of Fig. |2] indicate the strength 
of the break. To characterize the break strength we per- 
form a two-sided KS test on the best-fit SPL and BPL den- 
sity models. Weak breaks will have a high probability of be- 
ing drawn from the same distribution as the SPL profile (i.e. 
^KS ^ 1)^ while strong breaks will have a low probability 
(i.e. Pks ^ 0). One can see that the points deviating the 
most from a one-to-one relation between apocentre and break 
radius tend to have the weakest breaks. 

The bottom panel of Fig. |2]shows the average time of (stel- 
lar) stripping for each satellite (Tub, in Gyr) against the av- 
erage apocentre at z = 0. Satellites accreted at early times 
have smaller apocentres as the physical size (and mass) of the 
parent galaxy is smaller. Therefore, breaks give an indication 
of when a satellite was accreted and/or disrupted. However, 
a further complication is dynamical friction, as more massive 
satellites can sink to the center of the parent halo relatively 
quickly. The colors in the bottom panel of Fig. |2]indicate the 
satellite (stellar) mass. More massive satellites tend to have 
smaller apocentres (at stripping) even if they are accreted rel- 
atively recently. 

In Fig. [3] we show the radial velocity structure of three 
accreted satellites in HaloOS. The stripped material from the 
satellite in the left-hand panel shows an obvious break in the 
density profile which coincides with the apocentres of its star 
particle orbits. The middle panel shows a (massive) satellite 
accreted a long time ago (^ 9 Gyr ago) which is now well 
mixed in phase-space an d can be described by a single power- 
law density profile (cf. iJohnston et al.l l2008l) . Finally, the 
right-hand panel illustrates a relatively recent accretion rem- 
nant for which the density profile is poorly defined. In gen- 
eral, satellites accreted between 4.5-9 Gyr ago follow a BPL 
density profile. Satellites accreted a long time ago Tub si 9 
Gyr, can be well mixed in phase-space, and recent accretion 
events (Tub ^ 4.5 Gyr ago) are unrelaxed and their density 
profiles are poorly defined. 

We have found that the satellite apocentre sets the scale of 
the break radius, but what sets the 'strength' of the break? In 
the top panel of Fig. |4]we show the break strength (defined by 
a two-sided KS test between SPL and BPL profiles) against 
the spread of apocentres for stars belonging to an individual 
satellite. We find that stronger breaks have a smaller spread 
in particle apocentres. The spread in apocentres is related to 
the energy distribution of the star particles; the bottom left- 
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Figure 2. Top panel: The break radius of the stellar density profile 
against the average apocentre of the star particle orbits. The circle 
symbols indicate individual accreted satellite remnants from all 1 1 
stellar haloes, and the star symbols indicate their global values. The 
colors indicate the probability that the best-fit single power-law den- 
sity profile is drawn from the same distribution as the best-fit broken 
power-law density profile (i.e. using a two-sided KS test). The red 
colors indicate profiles which do not have a strong break in the stel- 
lar halo density. The dotted line shows a one-to-one relation. Bot- 
tom panel: The average time of stripping for each accreted satellite 
against the apocentre of the orbits. Satellites accreted at early times 
have smaller apocentres as the physical size of the parent halo is 
smaller. The colors indicate the satellite (stellar) mass. More mas- 
sive satellites are more affected by dynamical friction, so can sink to 
the center of their parent halo more quickly than less massive satel- 
lites. 

hand panel of Fig. |4] shows that the dispersion in energy of 
an accreted remnant is strongly correlated with the spread in 
particle apocentres. Note that these relations also hold for 
the global stellar halo properties (shown by the blue star sym- 
bols). Here, the spread in apocentre is calculated from all 15 
satellite accretion remnants. 

The bottom right-hand panel of Fig. |4] shows that more 
massive satellites ha ve a larger sprea d in energy (as expected 
analytically, see e.g. lJohnstonlll998l) . The colors indicate the 
average time of stripping for each satellite. Satellites accreted 
a long time ago (Tub ^ 9 Gyr) can have a large spread in 
energy, regardless of the satellite mass (cf. phase-mixed ex- 
ample in middle-panel of Fig. O. 

5. IMPLICATIONS FOR GLOBAL STELLAR HALO PROPERTIES 

In the previous section, we found that the density profile of 
stars stripped from an individual satellite is strongly related 
to the orbital structure of the stars. The average apocentre of 
the star particle orbits corresponds to the break radius of the 
best-fit BPL density profile. Furthermore, the strength of the 
break depends on the spread in energy (and hence apocentres) 
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Figure 3. Examples of individual satellite contributions to the global stellar halo. Three satellites accreted/disrupted at different times in Halo05 
are shown. The top panels show their radial velocities as a function of radius and the bottom panels show their density profiles. The dashed- 
line is the best-fit broken power-law model. In the left-hand panel, there is an obvious 'break' in the stellar density which coincides with the 
average apocentre of the star particle orbits. Satellites disrupted a long time ago (Tub ^ 9 Gyr) are often well phase-mixed and follow a single 
power-law profile. Recent accretion events (Tub < 4.5 Gyr) are often unrelaxed and have poorly defined (i.e. non-smooth) density profiles. 
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Figure 4. Top panel: The 'strength' of the broken power-law density model against the dispersion of star particle apocentres. The strength of the 
break is characterized by a KS test between single and broken power-law density models. A high probability that the best-fit single power-law 
model is drawn from the same distribution as the best-fit broken power-law model indicates a weak, if not non-existent, break in the stellar 
density profile. The black circles are individual satellites and the blue stars indicate the global stellar halo properties. Bottom-left panel: The 
dispersion in energy of star particles against the dispersion in star particle apocentres. A large range of star particle energies leads to a large 
dispersion in their apocentres, and hence a weaker break in the stellar density profile. Bottom-right panel: The dispersion in energy of star 
particles against satellite mass. Higher mass satellites have a larger range of particle energies. The colors indicate the time of stripping. 
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Figure 5. Examples of three stellar haloes. The black squares show the break radii of the 15 most massive accretion events against disruption 
time (left panels) and strength of break (right panels). The sizes of the square symbols indicate the fractional contribution of each satellite 
to the overall stellar halo (larger squares are more massive satellites). The blue dotted line and blue triangles indicate the overall stellar halo 
properties. The red stars indicate recent accretion events (Tub ^ 4.5 Gyr) where the density profiles of the remnants are poorly defined. These 
are placed at rj, = 90 kpc for illustrative purposes. The dashed black line indicates the mass-weighted break radius. This is calculated by 
computing the weighted mean of all 15 satellites' break radii. Haloes 14 and 07 (top and middle panels) can both be described by a broken 
power-law stellar density profile (albeit Halol4 has a relatively weak break: Pks ~ 0.99 and 02 — «i = 1-5) . In the latter case the accreted 
satellites all have similar apocentres (and hence break-radii) whereas the former halo is dominated by one massive satellite. HaloOS (bottom 
panel) has no obvious break (Pks ~ 1 and 02 — Qi = —0.2) and its accreted satellites have a wide range of apocentres. 
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Figure 6. Radial velocity dispersion profile for two example haloes. 
The black lines show the overall profiles, the red/blue lines show 
the profiles for star particles with metallicity greater thanAess than 
[Fe/H] = —0.9. The dashed lines show the profiles for the stars 
belonging to the most massive satellite accreted by the galaxy. The 
dotted lines indicate the break radii of each halo. The stellar halo of 
Halol4 is dominated by one massive satellite whilst HaloO? is built 
up from several similar mass satellites with similar apocentres. The 
metal-rich(er) velocity dispersion profile for Halol4 shows a more 
pronounced 'dip' near the break radius. There is little difference 
between the metal-rich(er) and metal-poor(er) profiles for HaloO?, 
indicating that the break is not dominated by a single, massive accre- 
tion event. 

of the stars. In general, more massive satellites have a larger 
energy spread. However, very early accretion events (Tub 9 
Gyr ago) which are well mixed in phase-space can also have a 
large spread in energy, regardless of the initial satellite mass. 

We now relate these findings to the global stellar halo prop- 
erties. The global stellar halo is a superposition of individual 
satellite contributions. Therefore, the overall stellar density 
break radius reflects the individual satellite contributions. For 
example, a strong break ensues when the accreted satellites 
all have similar apocentres (e.g. Halo07, middle panel of Fig. 
|5]). The halo can also be dominated by one massive satellite 
(e.g. Halol4, top panel of Fig. |5]l, and the overall break ra- 
dius reflects the apocentre of this massive satellite. In Fig. |5] 
we show the best-fitting break radius for individual satellite 
contributions versus average time of stripping (left-hand col- 
umn) and strength of break (right-hand column). The size of 
the black squares illustrates the satellite mass (larger symbols 



for more massive satellites). The blue triangles indicate the 
global parent stellar halo properties and the blue dotted line 
indicates the parent stellar halo break radius. We also show 
the 'mass-weighted' global break radius by the black-dashed 
line. This is computed from the individual satellite break radii 
weighted by their mass. The top two rows illustrate the two 
different scenarios mentioned above. In Halo 14 the overall 
stellar halo density profile is dominated by one massive satel- 
lite, whilst in HaloO? several satellites (of comparable) mass 
have very similar break radii. 

Are we able to discern between these two scenarios? Fig. 
|6] shows the radial velocity dispersion profiles for these two 
haloes. The black line is the overall profile, and the red and 
blue colors illustrate the profiles for metal-rich(er) and metal- 
poor(er) stars (arbitrarily split at the average metallicity of the 
stellar haloes [(Fc/H)] = —0.9). The profile for the most 
massive satellite contributing to the stellar halo (between 10- 
100 kpc) is shown by the black dashed line. In both HaloO? 
and Halo 14 there are 'dips' in the velocity dispersion profile 
near the break radius (indicated by the vertical dotted line). 
This is a result of the low radial velocities of the star particles 
at the apocentres of their orbits. However, in Halo 14 this dip 
is much more pronounced in the higher-metallicity stars. This 
is because a massive satellite, which is therefore relatively 
metal-rich, dominates the break in Halo 14. One can see that 
the profile for this massive satellite (dashed black line) closely 
follows the profile for the most metal-rich stars in the halo. 
There is little difference between the metal-rich and metal- 
poor material in HaloO?. This is because all of the satellites 
in HaloO? (low and high mass) have very similar apocentres, 
and there is no massive, dominating accretion remnant. 

Finally, we consider the case where there is no obvious 
break in the stellar density profile. In this case, the accretion 
events may have a wide range of apocentres (e.g. HaloOS, 
bottom row of Fig. |5]). As the apocentre is related to the time 
of accretion, this suggests that the accretion history has been 
prolonged. The bottom row of Fig. |5] shows that there have 
been significant satellite disruption events over a wide time 
interval: from ^11 Gyr ago to recent events only ~ 1 Gyr 
ago. In comparison. Halo 14 and HaloO? show little activity in 
the last 6 Gyr. 

In summary, the stellar halo density profile, and in partic- 
ular the presence or absence of a break in the stellar density, 
can give us an important insight into the accretion history of 
a galaxy. 

5.1. The Milky Way andM31 

In the last few years, several groups have found that the 
Milky Way stellar halo has a 'b r oken' density profile (e.g. 
iBell et alJl2008l:lSesar et al. l l2Mll:lDeasoirerai1l20TTbl) . with 
riy = 20 — 30 kpc an d ctin = 2.3 — 2.8 and aout = 
3.8 — 4.6. Furthermore, iDeason et al.l (1201 lal) found a dip 
in the line-of-sighQ velocity dispersion profile of the metal- 
rich(er) {\¥e/Y{\ > —2) blue horizontal branch stars (see Fig- 
ure 6 in lDeason et alj|201 lah . This dip in velocity dispersion 
occurs between 20-30 kpc, in the approximate radial range 
where there is a break in the stellar density profile. 

Our conclusions in the previous sections suggest that these 
observations may have several important implications: 

• The relatively close-by break radius of the Milky Way 

In the radial range of the IDeason et al.l )2011aD study, 10-50 kpc, the 
line-of-sight velocity is a good approximation for the radial velocity 
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(rh = 20 — 30 kpc) is likely caused by accreted satellites 
with small apocentres (at the time of stripping). 

• The cold velocity dispersion profile near the break ra- 
dius provides further evidence for the break-apocentre 
connection. 

• A metal-rich bias in the velocity dispersion profile sug- 
gests that the break in the Milky Way stellar halo could 
be dominated by a (relatively) massive satellite(s). 

• The density profile has only been measured out to 
50 kpc, so we cannot rule out more recent accre- 
tion events with larger apocentres. However, the rel- 
atively smooth stella r density profile out to 50 kpc 
dPeason et alJl20lTbh suggests early accretion events 
and then somewhat quiescent evolution. 

Recently, [Gilbert et al.l (120121) measured the density profile 
of the M3 1 stellar halo ou t to 100 kpc (see al s o earlier work by 
iGuhathakurta et al."2005', llrwin et alj|200l llbata et al.ll2007l 
and Courteau et al. 2011). These authors find that the stel- 
lar distribution can be described by a SPL with a ^ 3.3, 
and there is no evidence for a break in the density profile. 
As we found in the previous section, no obvious break in 
the global stellar halo profile indicates that the stars have 
a wide range of apocentres, which suggests that M31 has 
had a much more prolonged accretion history than the MiUcy 
Way (see e.g. HaloOS). This deduction, that M31 appears 
to have experienced a much more active accretion history 
than the Milky Way, is in good agreement with other in- 
dependent observations of the two galaxies: R elative to the 
Milky Way, M31 has a more disturbed disc (e.g. Brown et al.' 
2006h. a larger bu lge (e.g. [Pritchet & van den Bergh 1994; 



Durrell et al.l | 2004l). a younger and more me tal-rich halo pop- 



ulation (e.g. llrwin et ani2005t iKalirai et al.ll2b06i) . and more 
numerous tidal streams and surviving sa telhte galaxies (e.g. 
iKoch et al.ll2008tlRichardson et alJl20lTl) . 

6. CONCLUSIONS 

The stellar haloes of the Milky Way and M31 have con- 
flicting density profiles: the Milky Way has a broken pro- 
file, whereby the stellar density falls off more rapidly be- 
yond a break radius (rf, ^ 25 kpc), whereas, the stellar 
halo of M31 shows no obvious break and can be described 
by a single power-law. In light of these recent observations, 
we have studied the density profiles of stellar haloes - built 
solely by the accret i on of dwarf galaxies - drawn from the 
iBuUock & JohnstonI (I2005h suite of simulations, with the aim 
of understanding the contrasting stellar haloes of our local 
galaxies. 

We summarize our conclusions as follows: 

(1) The simulated haloes often have 'broken' stellar halo pro- 
files, where the density falls off more quickly beyond the 
break radius. However, some haloes do not have an obvi- 
ous break, and their density distribution can be described by a 
single power-law (e.g the MiUcy Way versus M31). 

(2) In the BJ05 simulations, the global stellar halo is a su- 
perposition of accretion products from the ^ 15 massive ac- 
cretion events. The density profiles of the stars that once be- 
longed to an individual satellite follow a broken profile, where 
the break radius corresponds to the average apocentre of the 



stars. However, material belonging to recently stripped satel- 
lites (Tub ^ 4.5 Gyr) have ill-defined density profiles, and 
very ancient accretion events (Tub ?J 9 Gyr) are often well- 
mixed in phase-space. 

(3) The location of the break in the stellar density is linked 
to the time of accretion/stripping and the mass of the satel- 
lite. Satellites accreted at early times have smaller apocentres 
as the physical size of the parent halo is smaller More mas- 
sive satellites can have smaller apocentres (and hence break 
radii) at the time of stripping, as the satellite can spiral into 
the center of the galaxy via dynamical friction. 

(4) The strength of the break in the stellar density depends on 
the spread of apocentres of the stars. For the global stellar 
halo, the break strength depends on the range of apocentres of 
its accreted components. Individual satellite accretion rem- 
nants can also have varying break strengths. More massive 
satellites have a larger spread in energy (and hence apocen- 
tres). Also, stars belonging to satellites accreted a long time 
ago can have a wide spread in energy as they are well mixed 
in phase-space today. 

(5) The global stellar halo density profile depends on the ac- 
cretion history of the galaxy. An obvious break in the overall 
density profile suggests that: (a) the accreted satellites have 
similar apocentres; or (b) one massive satellite dominates the 
stellar density in the applicable radial range. These two sce- 
narios could be distinguished observationally. The radial ve- 
locities of stars near apocentre are very low, thus the radial 
velocity dispersion is also low (i.e. shell-type structures). In 
case (b) above, the apocentres of stars that once belonged to a 
massive, and hence metal-rich, satellite dominate the break. 
Thus, a 'dip' in the velocity dispersion profile, the signa- 
ture of a build-up of apocentres, will be more pronounced in 
the metal-rich(er) material. Conversely, if the accreted satel- 
lites all have similar masses (and similar apocentres), then we 
would expect no metallicity bias. 

(6) The profiles of some haloes show no obvious break in the 
stellar density. Often, these haloes have a prolonged accretion 
history whereby satellites are accreted over a wide range of 
timescales (~ — 10 Gyr ago). Thus, the average apocentres 
of the accreted satellites, whose material now makes up the 
stellar halo, also have a wide range of values. 

Our investigation into the simulated BJ05 stellar haloes has 
provided some important insights into the accretion histories 
of the Milky Way and M31. The absence of a break in the 
M3 1 stellar halo suggests that this galaxy has had a prolonged 
accretion history, where the accreted satellites had a wide 
range of apocentres. However, the strong break in the Milky 
Way stellar halo (at = 20 — 30 kpc), in addition to the pres- 
ence of a shell-type feature in the relatively metal-rich BHB 
stars, suggests that the stellar halo break is dominated by the 
apocentre of a (relatively) massive satellite. 
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APPENDIX 

RELATING SATELLITE PROPERTIES AT ACCRETION TO Z = HALO STARS 

In this section, we verify some of our deductions from the star particle properties at z = 0, by tracing back the accreted 
satellites whose stripped stellar material make up the stellar haloes today. 

We follow the orbits of satellites from Halo02 and Halo07. In the left-hand panels of Fig. [7] we show three examples of 
satellites from Halo07. These are the three most massive contributors to the stellar halo within 100 kpc of this galaxy. In the left 
columns we show the evolution of radius with time for the satellites, and the right columns show the evolution of total satellite 
mass with time. The red dashed line indicates the average time when stars become unbound (< Tub >) from the satellite. The 
dotted black line indicates the average apocentre of the star particles computed in Section|4] and the solid blue line indicates the 
approximate break radius in the density profile of the stripped stellar material. This figure shows that the break radius we measure 
at z = from the stars is indeed coincident with the apocentre of the host satellite galaxy near the time of stripping. In Section 
m we deduced this by estimating the apocentre of star particle orbits using Equation[T] 

In the right-hand panels of Fig|2]we relate some of the initial satellite properties to those we infer from their stripped stars 
at z = 0. The top-left panel of this figure shows the relation between accretion time (of the satellite) and the average time at 
which the star particles become unbound. The stars are stripped a few Gyr after satellite accretion, but there remains a positive 
correlation between these two timescales. The bottom-left panel shows that the difference between these timescales is related to 
the mass of the satellite galaxy. Due to dynamical friction effects, more massive satellites spend a shorter amount of time orbiting 
the galaxy before the stellar material is stripped. In the top-right panel we show satellite accretion time against initial satellite 
apocentre (red asterisks) and average star particle apocentres at z = (filled black circles). In Section ID we noted that larger 
apocentres correspond to more recent accretion events as the physical size, and mass of the parent galaxy is larger, and hence 
more distant satellites can be captured. The strong trend between initial apocentre and accretion time reinforces this statement. 
Although the trend is slightly weaker, we see that the apocentres of the stripped stars at z = are still an indication of when 
their host satellite was accreted. The final apocentre of the satellite (and hence stripped stars) is generally reduced from the initial 
apocentre via dynamical friction effects. In the bottom-right panel of the figure we show that the difference in apocentre radius 
is strongly related to the satellite mass. 

In summary, we find that our deductions made from the z — halo stars in the main section of the text, are verified when we 
trace the accreted satellite properties back in time. 
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Figure?. Left hand panels: Three example satellite orbits from HaloO?. The radius (left columns) and mass (right columns) are shown as a 
function of time. The satellites are followed from their initial apocentre (at accretion) until the satellite is completely disrupted. The dashed red 
line shows the average time at which the stars belonging to the satellite are stripped. The dotted black lines show the average apocentre of the 
star particles (at stripping) calculated in Section|4] The solid, blue line shows the estimated break radius of the stripped stellar material. Here, 
we confirm our deductions in the main text; the break radius is coincident with the apocentre of the satellite near stripping. Right hand panels: 
The top-left panel shows that the average time at which stars are stripped from the satellite is correlated with the time at which the satellite is 
accreted. The top-right panel shows that the accretion time of a satellite is correlated with both the initial apocentre (red asterisks), and the 
final apocentre (filled black circles). The bottom two panels show the effect of dynamical friction. We show the orbital time (left panel) and 
reduction in apocentre (right panel) against satellite mass. More massive satellites are stripped quicker, and have a larger difference between 
initial and final apocentre. 



